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Abstract 

Circularly rotating axisymmetric perfect fluid space-times are in- 
vestigated to second order in the small angular velocity. The condi- 
tions of various special Petrov types are solved in a comoving tetrad 
formalism. A number of theorems are stated on the possible Petrov 
types of various fluid models. It is shown that Petrov type II solu- 
tions must reduce to the de Sitter spacetime in the static limit. Two 
space-times with a physically satisfactory energy-momentum tensor 
are investigated in detail. For the rotating incompressible fluid, it is 
proven that the Petrov type cannot be D. The equation of the rota- 
tion function uj can be solved for the Tolman type IV fluid in terms of 
quadratures. It is also shown that the rotating version of the Tolman 
IV space-time cannot be Petrov type D. 



1 Introduction 

Following Schwarzschild's discovery of the static incompressible interior so- 
lution in 1916, relentless efforts have been made to find a rotating general- 
ization. Improved equations of state for the perfect fluid were imported in 
general relativity, leaving, however little clue for how to achieve this goal. 
Less than that, the supply of even the physically acceptable nonrotating so- 
lutions is scarce |l[]. The common approach to solve the field equations of 
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the problem is to make some a priory assumption on the properties of the 
desired spacetime, i. e., make an ansatz, and hope that solutions with the cor- 
responding property do exist. More often than not, after considerable effort 
one gets a highly unphysical solution in this procedure, if any solution turns 
out to exist at all. A natural way to decide the acceptability of the various 
ansatze is to check whether or not they remain valid in the slowly rotating 
limit. Since static spherically symmetric spacetimes are either Petrov type D 
or 0, a widely studied ansatz is to assume that the rotating fluid also belongs 
to some special Petrov class. The aim of this paper is to investigate whether 
or not these classes contain physically acceptable perfect fluid spacetimes. 
It is natural to assume, that such configurations have a well behaving slow 
rotation limit. 

We investigate axisymmetric stationary perfect fluid spacetimes in cir- 
cular rotation, i.e., with the fluid velocity vector lying in the plane of the 
timelike and angular Killing vectors, d/dt and d/dip respectively. The met- 
ric of the spacetime is written in the form 



d^^ + sin^ ^{dip - Ludt)^ , (1) 



where X, Y, Z and u are functions of the coordinates r and We choose 
the time-translation Killing vector d/dt such that it becomes asymptotically 
nonrotating at spacelike infinity. With this choice, the angular velocity fl 
of the fluid is defined by the components of the fluid velocity vector as 
follows, 

= nu' . (2) 

We further assume that the fluid is in rigid rotation, i. e. is a constant. This 
Q is the parameter which is small in the slow rotation limit, and following 
Hartle's work ||^, we will expand quantities in powers of Q. For rigidly 
rotating space-times, a comoving coordinate system, where u'^ = 0, can be 
arranged by a linear transformation of the angular coordinate if — > if — fit. 
Since in this work we are interested only in the interior fluid region, we will 
perform our calculations in the comoving system. 

We focus on the Petrov types that a slowly rotating fluid ball might 
have. We use a tetrad, constructed specifically for this task, representing the 
space-time to the desired order in the angular velocity Q. In this framework, 
definite statements can be made about the Petrov type of the field of a given 
order in Q, or about the impossibility of certain Petrov types to any order 
in Q. Nevertheless, the slow rotation approximation can be applied to give 
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information about numerous physical properties, even about the possible 
existence of an asymptotically flat exterior vacuum region [Q. 

Since the system behaves the same way under a reversal in the direction 
of rotation as under a reversal in the direction of time, when expanding the 
components of the metric (|l|) in powers of the angular velocity parameter 
Q one finds that uj contains only odd powers, while X, Y and Z contain 
solely even powers of fl. In this paper we are interested in effects of up to 
order, and hence we consider the function u to be proportional to the angular 
velocity parameter Q, while we allow Q"^ terms in the diagonal components 
of the metric (|I]). The metric of a slowly rotating fluid ball can be written 
in the form [Q] 

ds^ = {l + 2h)X^dt'^ - {l + 2m)Y^dr^ 

-{l + 2k)Z^[d^'^ + 8111^^ {dip- iudtf] . (3) 

Here the functions X, Y and Z depend only on the radial coordinate r, 
determining the spherically symmetric basis solution, while u, h, m and k 
are functions of both r and i}. The potential lj is small to first order in the 
angular velocity fl, but h, m and k are second order small quantities. There 
are two minor differences between our metric form and the corresponding 
formula in the paper of Hartle [@. The first is that we do not use the 
radial gauge Z = r, since that choice is technically inconvenient for certain 
spherically symmetric perfect fluid exact solutions. The second difference is 
that our definition of the second order small quantity m differs by a zeroth 
order factor. This choice is only to make our equations shorter, especially 
because the local mass function in the denominator of Hartle's definition 
takes a more complicated form in a general radial gauge Z ^ r. 

Let us consider a freely falling observer, with velocity vector v^, who 
has zero impact parameter and consequently zero angular momentum. In 
the coordinate system where d/dt corresponds to the timelike Killing vector 
which is nonrotating at spacelike infinity, v^p = and the function oj agrees 
with the angular velocity v"^ /v^ = g'^^/g^^ = —g^t/g<fiifi of this observer. Hence 
uj can be interpreted as the rate of rotation of the local inertial frame with 
respect to the distant stars[0, or in other words, the cumulative dragging of 
initial framesQ]. In comoving coordinates, the rotation potential u represents 
the coordinate angular velocity of the fluid element at (r, measured by a 
freely falling observer to first order in Q. 

We next briefly recapitulate Hartle's arguments leading to an unambigu- 
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ous choice of the coordinate system. For the full details, cf. [Q. We start 
from a known spherically symmetric perfect fluid solution, which is described 
by the functions X, Y and Z. Calculating the field equations to first order 
in the small angular velocity parameter Q, we get only one independent re- 
lation P The {t, if) component of Einstein equation gives a second order 
partial differential equation for uj. In the comoving coordinate system this 
equation takes the form 



X d ( duj\ XZ 
+ 4- 



drXXYdr Y 



d I dZ 



1 ^duj\ ^ 

"+;h?^M r°- 



dr \XY dr 

' (4) 

Expansion of @) in vector spherical harmonics yields the angular behavior 
of the solution in the form 



id = ^UJi 

1=1 



1 dPi{cos^) 
sin dd 



(5) 



where the functions uji depend only on the radial coordinate r, and Pi is the 
Legendre polynomial of order I. The equations for the coefficients ui with 
different values of I decouple. Taking into account the matching conditions 
at the surface of the fluid ball, one can show that for / > 1 the functions uji 
cannot be regular both at the center of the fluid and at infinity. It follows 
from the asymptotic fiatness of the exterior spacetime region that uj cannot 
depend on the angular coordinate d. Consequently, the rotation potential 
a; is a function of the radial coordinate r alone, even in the fluid region, 
satisfying a second-order ordinary linear differential equation. 

Since the expansion of uo in the angular velocity parameter cannot contain 
terms, the solution of (Q) will remain valid to second order as well. After 
solving the first order condition for the function uj, proceeding to second 
order in the angular velocity parameter f2, the components of the Einstein 
equation give a coupled linear inhomogeneous system of partial differential 
equations for the functions m and k. The inhomogeneous terms in the 
equations are proportional to uj"^ and its derivatives. The solution of this 
system can be written in the form of an expansion in spherical harmonics, 

00 

h = J2hiPi{cos'd) , (6) 

and similarly for m and k. However, quantities with different / decouple, and 
the equations for I > 2, being homogeneous, do not include u. As a result. 
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all hi, mi and ki must vanish for Z > 2, since otherwise they would correspond 
to a static but not spherically symmetric configuration in the ui = case. 
Thus the second order small functions can be written in the form 



m 



h 



k 



ho + h2P2{cos{}) , 
mo + m2P2(cos'i9) , 

k2P2{cOS^) , 



(7) 



where ho, h2, mo, m2 and /c2 are functions of r. The freedom in the choice 
of radial coordinate was used to set the independent part of k to zero. 

In Ref. 0, Hartle writes out the detailed form of the equations describing 
the second order rotational perturbations in the Z = r gauge, and gives a 
procedure to determine the binding energy, the baryon number change, and 
the ellipticity of the fiuid surface. However, this procedure involves numerical 
integration of a system of ordinary differential equations. Instead, in our 
paper, we focus on what can be said about the general physical properties of 
slowly rotating fluid bodies by analytical methods, without trying to solve 
the perturbation equations. There is an extensive literature on numerical 
simulations of slowly rotating bodies for various types of fluids. These include 
neutron stars and supermassive stars , incompressible fluids ^ , polytropes 
and realistic neutron matter equations of state 0. 

In Sec. 2., we establish a tetrad formalism for the Petrov classiflcation. 
We also prove a theorem showing that physically realistic slowly rotating 
perfect fluid spacetimes cannot be Petrov type II. Case studies for various 
equations of state are presented in the rest of the paper. In Sec. 3., we 
establish a theorem on incompressible fluids. This theorem indicates that 
circularly rotating states should be found in the algebraically general class. 
For Tolman IV fluids, in Sec. 4., we flnd the rotation function u in terms of 
quadratures, and we show that these fluids cannot be Petrov type D. 

The rotation function u can be written down in terms of elementary 
functions for the rotating Whittaker space-time. Although the equation of 
state of the Whittaker fluid implies that the density decreases towards the 
center of the fluid ball, this class merits special attention. Among the rotating 
states of the Whittaker fluid, there is the exactly known Petrov type D 
Wahlquist solution. The rotating Whittaker fluid is dealt with in and 
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2 The Petrov types 



To quadratic order in the angular velocity parameter f2, the nonvanishing 
components of a comoving tetrad can be chosen for the metric (|^) as 

ei = ^^, 62=^—, 63 = ^— smT?, (8) 

The correctness of these expressions to the required order can be shown by 
checking that (e'^e^g^u) = diag(l, —1, —1, —1) up to fl^ terms, where Roman 
and Greek labels denote tetrad and spacetime indices, respectively, and g^^, 
are the components of the metric (|^). 

Using the tetrad components (|), we will compute the Ricci rotation co- 
efficients, the tetrad components of the Riemann, Weyl and Einstein tensors, 
and also some expressions, which are vanishing for certain Petrov types. All 
these quantities are polynomial expressions in the metric components gfj^^, 
the tetrad vector components e^, and their coordinate derivatives. Since the 
tetrad vector components will be multiplied only by terms which have a reg- 
ular f2 = limit, we get the correct results to second order in Q even if we 
use expressions for the tetrad, which are only correct to the same order.Q 

The electric and magnetic curvature components are defined in terms of 



the tetrad components of the Weyl tensor as follows |]TT 



El — ClOlO , E2 — C2O2O ; — C102O ; 

-ni — L/ioiO , -"2 — <-/2020 ; -"3 — tvl020 



^An alternative approach would be to consider temporarily (|3|) to be valid to arbitrary 
order in the angular velocity f2, calculate the exact tetrad vector components, the Riemann 
tensor and other necessary quantities, and finally expand the result to second order in fl. 
This procedure gives exactly the same results as the simpler one, but because of the 
appearance of complicated denominator terms it would be a tedious work even for modern 
computer algebra systems. An even more cautious way of calculation, leading to the same 
results once again, would be to compute the coordinate components of the Riemann tensor, 
and take the tetrad components only in the end. However, this approach would have a 
similar problem with the components of the inverse metric g^"^ , namely, whether or not 
one calculates them only to second order from the beginning. 
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Because of the symmetries of the configuration, these six components are the 
only independent components of the Weyl tensor. A Newman- Penrose frame 
is readily defined by I = (cq + 63) / a/2, n = (cq — 63) / v^, m = (ei + 162) / \/2. 
The nonvanishing Weyl spinor components are 

^0 = ^ [E2 -El- 2H3 + i{Hi-H2- 2^3)] , (9) 

^2 = ]^[Ei + E2-i{Hi + H2)] , (10) 

^4 = \[E2-Ei + 2H^ + i{Hi-H2 + 2E^)]. (11) 

The Petrov type of the spacetime can be determined by studying the prop- 
erties of the Weyl spinor components. If \E'o = ^2 = ^4 = the spacetime 
is conformally fiat, consequently it is the nonrotating interior Schwarzschild 
solution|T^. If \E'o = ^4 = but ^2 7^ 0, one can apply a tetrad rotation 
around the vector 1, which makes m into m + al. Then the new nonzero 
Weyl spinor components become \E'4 = 6(a*)^\E'2, ^'3 = 3a*\I'2 and \E'2 = '^2-, 
and the Petrov type is determined by the number of the distinct roots of the 
equation = ^46'' + 4^36^ + 6^2^^ = 6\[^2&^(a*& + 1)^- This equation has 
two double roots, which shows that the spacetime is of type D. If \E'4 7^ the 
Petrov type is determined by the number of the distinct roots of the algebraic 
equation 

^4&^ + 6*2fe^ + ^0 = (12) 

for the parameter h. If \E'4 = but \E'o 7^ one can apply a reversal in 
the 63 direction to interchange the frame vectors 1 and n, and consequently 
interchange and \E'o as well. Type III solutions are excluded because ^3 
and \&i vanish. The type is N if = ^2 = but \E'4 7^ 0. It has been 
shown in |ll| that all axistationary type N perfect fluid solutions can be 
interpreted as vacuum solutions with a negative cosmological constant. In 
the case \E'4 7^ the Petrov type is II if and only if \E'o = but \E'2 7^ 0. 
Assuming ^^4 7^ 0, the Petrov type is D if and only if 9\E'2 = ^0^4 7^ 0. For 
our generic comoving tetrad, the real and imaginary parts of the equation 
9^2 = ^0^4 yield 

2El + 5E1E2 + 2^2^ - Ej - 2Hl - 5H1H2 - 2Hl + Hi = , (13) 
AEiHi + 5E1H2 + 5E2H1 + AE2H2 - 2^31/3 = . (14) 

These conditions characterize slowly rotating fields. For finite angular veloc- 
ities, the conditions for higher-order terms may not be met any more. By the 
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generic properties of the power series expansion of different iable functions, 
however, any statement about the nonexistence of a certain type slowly ro- 
tating fluid will hold in the exact sense. 

In the spherically symmetric limit, when = 0, it is easy to check that 
\l'o = ^1/4 = — 3\l/2, independently of the equation of state. There are only two 
possible Petrov classes then. If all are zero, the metric is the conformally 
fiat Petrov type incompressible interior Schwarzschild solution. All other 
spherically symmetric perfect fluid spacetimes are Petrov type D, since 9\l/2 = 
\l/o^4 7^ holds for them. This has important consequences for rotating 
Petrov N and Petrov II fluids, and also for the \l/o = ^4 = Petrov D 
subclass. Since \l/4 or \l/o vanishes in all of these cases, in the i7 = limit 
all the Weyl spinor components must go to zero, and consequently these 
spacetimes must reduce to the incompressible interior Schwarzschild solution 
in the nonrotating limit. This is already a severe limitation, but as we will 
see shortly, higher order conditions on the slowly rotating fluid state will 
pose further restrictions, which make these classes irrelevant for the study of 
rotating isolated bodies. 

Lemma 1: Circularly and rigidly rotating perfect fluids with \l/o = 
must reduce to the de Sitter space-time in the slow-rotation limit. 

Proof: Using @, the real and imaginary parts of the condition \E'o = 
give 

^2 - El - 2/^3 = , (15) 
H1-H2-2ES = . (16) 

We can prove the lemma without assuming the existence of an asymptoti- 
cally flat region, i. e. when u may also depend on We denote the tetrad 
components of the Einstein tensor by Gab- As we have discussed in the pre- 
vious paragraph, the spacetime must reduce to the incompressible interior 
Schwarzschild solution in the nonrotating limit, and hence the metric is of 
the form (|^) with 

X = A — cosr , 

Y = R, (17) 
Z = Rsinr , 

where A and R are constants, satisfying R > and 1 < A < 3. It is also 
possible to get this result directly by solving the zeroth order parts of the 
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equations (|T3|) and Gu = 6*22 • A linear combination of the first order parts 
of Go3 = and (0) yields a differential equation for u without ^ derivatives. 
The general solution of that is 



LU = 



( 



cos(r) + 1 



2A 



1 /i(^) + 



cos(r) — 1 



2A 



1 /2(^) , 



(18) 



where fi{^) and f2{'&) are some functions of the 'd coordinate. This uj is 
regular at the center of the fluid r = only if f2{'&) = or if A = 0. 
Supposing that A is nonzero, substituting into the leading cos^ r terms 
yield 



The solution of the second factor is proportional to 1/ sin^i), which is di- 
vergent at the rotation axis = 0. Consequently we must have A = 0. 
Substituting back again to (|T^) and (|TB|) we get that the derivative of fi{^) 
must be zero, and consequently is a constant. Since A = 0, this is the 
de Sitter space-time in a rotating coordinate system. The pressure p and 
density /x are constants, p = — /i = 3/R. Thus the equation of state violates 
the weak energy condition. 

The Petrov type is II in the following two cases: if and only if \E'o = 
but \l/2 and \l/4 are nonzero, or when \l/4 = but \I'2 and \I^o are nonzero. In 
the \l/4 = case we change the components of 63 to —1 times those in (|^) in 
order to exchange \l/o with \l/4. From our Lemma, it follows that 

Theorem 1: Circularly and rigidly rotating perfect fluids of Petrov type 
II must reduce to the de Sitter space-time in the slow-rotation limit. 

Our results also show that a slowly rotating perfect fluid with an equation 
of state satisfying the weak energy condition cannot be of Petrov type II. 

Lemma 1 also applies to the \l/4 = \l/o = but \l/2 7^ Petrov type D 
subcase. This shows that all physically acceptable rotating Petrov type D 
solutions must be in the Q^'l = ^0^4 7^ class. Since this case seems to be 
too complicated for general investigation, we proceed with studying perfect 
fluids with specific equations of states. 




(19) 
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3 Incompressible fluid 



The Schwarzschild metric of a non-rotating incompressible perfect fluid ball 
is described by (|17|) . The pressure and density are 



P 



3 cos r ~ A 



cosr 



(20) 



We assume the existence of an exterior asymptotically vacuum region, 
in which case u cannot depend on 6. Calculating to second order in the 
small angular velocity parameter Q, the {t, ip) component of the Einstein 
equations gives the condition 



(A — cosr) smr— — — (3 cos r — 4^4 cosr + 1 — AAusmr = 



(21) 



There are also three other second-order field equations involving the functions 
h, m and k. 

Using the tetrad (§), we get that the vorticity is linear in the angular 
velocity parameter f2, 



002 



COS r — A 
sin d 



2{A 



cosr 



2ijj\A cos r — 1) + sm r (A — cos r) — 

or 



(22) 
(23) 



The changes of all other rotational coefficients are small to second order. 

The electric part of the Weyl tensor is quadratic in the angular velocity 
parameter f2, and the magnetic part has only linear terms, 



cos?? 



R (cosr — A) dr 

sini} 
2R{A - cosr)2 
LUo COS r 

+ 



Ho = H^ 

I 2 



(24) 



dio 



cosr (cosr — A)— — h 2ylci;sinr 
dr 

ijj sin-i? 



R sin r R sin r (A — cos r) 



(25) 



In getting ifa, we used (|2l|) to eliminate the second derivative of uj. 

Since the other subcase is excluded by Lemma 1, the Petrov type can be D 
if and only if Eqs. ([13D and (14) hold. For the slowly rotating incompressible 
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fluid the only condition they give is iJ| — 2H'f — 2iJ| — 5H1H2 = 0, which 
imphes 

cosrfcosr — A)— — \-2ujAsmr = 0. (26) 
dr 

The general solution of this equation is 

/A \ 2 

u; = C^-l) , (27) 



Vcosr 

where C is a constant. However, this u is not a solution of (^1]) , and hence 
we have proven: 

Theorem 2: A slowly and circularly rotating incompressible perfect 
fluid spacetime with an asymptotically flat vacuum exterior cannot be Petrov 
type D. 

4 Tolman fluid 



A spherically symmetric perfect fluid space-time has been given in |]T3| . The 
metric is described by 

^2 ^ B^r' + A') 



2 + 2r^ 2 

^ = (i?2_^2)(^2+^2)^ ' (28) 

Z = r , 

where A, B and R are constants. The equation of state is quadratic in the 
pressure, 

R^{A^ + 2i?2)/i = 4i?MV + R^{2R^ + 13A^)p + 6(i?^ + 2^^) . (29) 

The metric has a pleasingly simple form and describes an isolated fluid body 
with a regular center, outwards decreasing density and pressure and sublu- 
minal sound speed for appropriately chosen parameters 

We assume the existence of an asymptotically flat exterior vacuum region, 
in which case the function uj does not depend on and the (t, Lp) component 
d^) of the Einstein equations gives 

(r^ - (2r2 + ^ + Ua^ - AA'— + 8r' - GrR'] ^ 

(XT' \ f J QjT' 

+4 {2R^ + A^) = . (30) 
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The general solution of this is 



^2r^ + 



C\ exp - 



+C2 exp 




(31) 



where C\ and C2 are constants and 



+A^ 
P3 



,10 



lOR' + 13A']r 



A* + AA^R^ - AR 



+ (6A^ + 4A^R^ - 8R' 
+ R^A^ [A^ - 6R^) - 2A^R'' (32) 



R\^ 
2 



R^ + A^] i^R^ + AA 



2r^ + A' 



(33) 



(r^ + 4r4i?2 ^ 2AV^ + r^A^ + 2r''A^R^ - A^R^) . (34) 

Since the \E'o = ^4 = subcase is excluded by Lemma 1, the rotating fluid 
can be Petrov type D only if equations (|T^) and (|l^) are satisfied. Using the 
tetrad (||) we get that (|1^) holds identically to second order in the angular 
velocity parameter VL. After substituting for the second derivative of uo from 
the Petrov D condition (|T3|) yields 

2 




+ Ar{r^-R^){A^ + 2R^)u 



dr 



(35) 



+4^2 V 



^2 ^ 6ij2^2 



2R^){A^ + 2r^ 



-{■1712 - /i2) = 



A2 + 2r2 A2(r2 + A2) 

Unlike in the incompressible case, this condition now involves the second 
order small quantities 7712 and /i2- The reason for the easier treatment of the 
interior Schwarzschild solution was that because of its conformal-fiat nature, 
both the electric and magnetic parts of the Weyl tensor were vanishing in the 
nonrotating limit. To be able to decide about the Petrov type, in general, one 
has to consider those parts of the Einstein equations, which are second order 
small in the angular velocity parameter Q. There are three such equations 
for the Tolman IV fiuid. The pressure isotropy condition 6*22 = ^33 takes 
the form 



yl2^4^^2_^2)(^2^2r2 



duj\ 
dr J 

+QB^R^{A^ + 2r^)\m2 + h 



) -4AV(A2 + 2i?2)a;2 
. 



(36) 



12 



The condition G12 = gives 



r(A^ + r^) + k2) - A^h2 - (A^ + 2r^)m2 = , (37) 
dr 

while the P2{cos'd) part of the another pressure-isotropy condition Gu = G22 
yields a more lengthy second order differential equation 



3BMr^-R^)iA^ + 2r' 



d\k2 + 112) / .2 , o 2\dm: 



dr"^ dr 



+35V(4r6 + Ah^ - 2R'r^ + 2R'A\^ + A^R^)^ (38) 

dr 

+3B^r\A^ + 2RyA^ + r'')^ 

dr 

-12B^R\A^ + 2ry{h2 + k2) + Ah\r^ - R^A^ + 2r^) = . 

In the following, we use equations (|35|) and ( p6D to eliminate /i2 and m2, while 
we employ ( PDD to express the second derivative of the rotation potential ui. 
Comparing Eq. (^8]) with the derivative of (^) we get another first order 
differential equation. Eliminating the first derivative of ^2 by (|37|) , we get an 
equation that can be solved algebraically for /c2- Substituting ^2 back to (|37D 
again, we obtain an equation containing only uo. Taking the r derivative of 
this equation, we readily find that it is not consistent with (0). Hence we 
have proven 

Theorem 3: A rotating perfect fluid spacetime that reduces to the Tol- 
man IV solution in the static limit and can be matched to an asymptotically 
flat vacuum exterior cannot be Petrov type D. 
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